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magnetic molecular complexes with low-energy singlet-triplet spin gap and show that the number 
of assisting phonons may be changed by varying the Zeeman splitting of excited triplet state. As 
a result, the structure of multiphonon Kondo resonances may be scanned by means of magnetic 
field tuning. 



INTRODUCTION 

Single electron tunneling through molecular bridges in nanodevices is inevitably accom- 
panied by excitation of vibrational modes. Vibration-assisted processes usually manifest 
themselves in tunneling through nanodevices as phonon satellites, which arise around main 
resonance peaks (see, e.g., [l|, y] and references therein). However, phonon assistance may 
induce resonance peak due to interplay with magnetic degrees of freedom in transition metal- 
organic complexes (TMOC). Appearance of Kondo-type zero bias anomaly in tunneling 
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through TMOC was predicted recently [3|. In this paper we consider multiphonon processes 
in Kondo tunneling in presence of magnetic field. We discuss the fine tuning effect of mag- 
netic field on the Kondo tunneling induced by multiphonon processes in a situation, where 
the ground state of TMOC with even occupation is a spin singlet. In this case the Kondo 
resonance in tunneling arises when the phonon emission compensates the gap between one 
of projections of the excited high-spin and the ground zero spin state of TMOC. Then the 
Kondo-like effect arises due to singlet /triplet transitions, which can be treated as effective 
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spin-flip processes |4| 



MODEL 

Phonon-assisted electron tunneling through TMOC as well as tunneling through other 
nanoobjects (quantum dots, nanotubes, ets) is described within a framework of Anderson 
model supplemented with the terms describing vibrational degrees of freedom and their 
interaction with electron subsystem 2|, y, |5j 

H = Hd + Hi + Hturi + Hyib + He-vib ■ (1) 

Here Hd stands for the electrons in the 3(i-shell of TMOC. It includes strong Coulomb and 
exchange interaction, which predetermine the spin quantum numbers of the corresponding 
electron configuration 3(i" {n is even)). Hi contains electrons in the metallic leads playing role 
of source (s) and drain (d) in the electric circuit, Htun is responsible for electron tunneling 
between TMOC and leads, H^b describes vibrational degrees of freedom in TMOC, and the 
interaction between electronic and vibrational subsystems is given by the last term ifg-mfe- 
Following the approach developed in [2, l3| we choose the phonon-assisted tunneling as a 
source of this interaction and represent vibrational subsystem by a single Einstein mode Q. 



Then three last terms in ([T]) are written as 

H.^b + Htun + ^e-«fe = hQh^h + J^ [te-"-^'' ''^ c[^J. + H.c.) . (2) 

aka 

Here the operators h stand for phonons, Ckau, d^ describe electrons in the leads a = s,d 
and the d-shell of TM ion, respectively, A is the electron-phonon coupling constant. When 
deriving ([2]) we assumed that the coupling is strong enough, so that multiphonon processes 



are treated exactly by means of the Lang-Firsov canonica' 
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transformation (cf. [2|). Another 
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canonical transformation of Schrieffer- Wolff (SW) type [7|] excludes -&<„„ from the Hamilto- 
nian ([1]) and maps it on an effective Hilbert subspace with fixed (even) number of electrons, 
singlet ground state and low-lying triplet excited state of TMOC The effective SW-like 
Hamiltonian is 

H,ff = - AS^ + hS, + Hi + JtS-s + JrR-s + H„b (3) 

(see [3|, l6|). Here Hd [the first two terms in ([3])] is represented only by spin degrees of 
freedom, A = Et — Es is the energy of singlet-triplet (S-T) transition, S, Sz is the spin 
operator and its 2;-projection, h = ^bqB is the parameter of Zeeman splitting in exter- 
nal magnetic field B. The electron spin operator is given by the conventional expansion 
s = I Ylikk' Tliaa' (^ka'^'^cr'Ck'a' where T is the Pauli vector. Unlike the conventional SW 
Hamiltonian, our Hamiltonian ([3]) contains one more vector R, which describes three com- 
ponents of S-T transitions [6]. This vector will be specified below. The electron-phonon 
interaction is now built into the effective exchange constants Jt = t^e"'^^'^^ "''^SEt and 
Jji = t'^e''^^^'' ^^^ /SEs- Here 5Et, 5Es are the energies of addition of an electron from the 
leads to the TM ion in a triplet and singlet states, respectively. 

MULTIPHONON PROCESSES IN KONDO COTUNNELING 

The essence of the mechanism of phonon-assisted Kondo cotunneling as it was formulated 
in Ref. 3|], is that the spin excitation energy gap A which quenches the conventional 
Kondo effect in a nanoobject with even occupation, may be compensated by the energy of 
virtual phonon emission/absorption processes. As a result the zero bias anomaly (ZBA) 
arises in tunnel conductance in spite of the zero spin ground state of the nanoparticle. This 
mechanism implies fine tuning | hfl — A | < E^ , where E^ is the characteristic Kondo energy, 
which is rather restricttive condition. To make situation more flexible, we address here to 
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FIG. 1: Single phonon con- 
nects singlet with spin 1 
projection of triplet. 



FIG. 2: n-phonon processes 
connect singlet with spin 1 
projection of triplet (nhQ < 

A). 



FIG. 3: (n + ?Ti)-phonon 
processes connect singlet 
with spin 1 projection of 
triplet {{n + l)hn > A). 



the case of strong electron-phonon interaction and apply magnetic field as an additional 
tuning instrument. The main ideas are illustrated by Figs. 1-3. The Zeeman term in the 
Hamiltonian ([3]) is responsible of splitting of the triplet state |T^) into 3 components with 
spin projections fi = 1, 0, 1. If the condition A — hQ — h < Ek is satisfied (Fig. [1]), then 
the states \S),\Ti) form effective vector operator R with components R^ = \^\Ti){S\, 
R'^ = \^\S){Ti\, Rz = |Ti)(Ti| — 15) (S'L which acts effectively as a spin 1/2 operator 



and enters the SW Hamiltonian 



see 



for further details). The term JrR. ■ s is 



responsible for Kondo-type resonance tunneling in accordance with the mechanism proposed 
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21 for zero magnetic field, where all three components of spin S = 1 are involved. If the 
resonance condition is fulfilled for n-phonon processes, A — h — nhQ < Ek, then the Kondo 
tunneling is assisted by virtual excitation of n-phonon "cloud" (Fig. [2]). If the condition 
A + h — {n + m)hQ < E^ is valid, then the opposite spin projection \Ti) is involved in 
Kondo tunneling, and the vector R has the components R~^ = ^/2\S){Ti\, R^ = v^|Ti)(5'|, 
Rz = \S){S\-\T,){T-,\. 

To find the contribution of phonon-assisted processes in Kondo tunneling, one should cal- 
culate the exchange vertex 7/1, which renormalizes the bare vertex J^ due to phonon emis- 
sion/absorption processes and logarithmically divergent parquet insertions. This dressed 
vertex is shown in Fig. HI To use the Feynman diagrammatic technique, the spin operators 
are represented via effective spin-fermion operators [f\. The wavy line corresponds to a 
single-phonon propagator in the case shown in Fig. [H Then the straightforward calculation 




FIG. 4: Phonon and parquet corrections to the vertex 7^. Solid and dashed hnes denote spin- 
fermion and conduction electron propagators, all parquet series are incorporated in the insertion 
shown by the square box. The multiphonon propagator is shown by the wavy line. 
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similar to that presented in [3] gives for the corresponding vertex 

D 



7hm 



pj^log 



maxlksT, h, | A — hQ — h\] 



l-{pJn)'\og' 



D 



(4) 



^niax[/cBT, /i, | A — hVt — h\\ 

Here T is the temperature, D is the effective width of the electron conduction band and p 
is the density of states on the Fermi level. The tunnel transparency T is proportional to 7^, 
and it is seen from (j4j), that the Kondo peak arises as a ZBA in T, provided | A — Ml — /i| ~ 



Ek ^ B. If the multiphonon processes are involved in accordance with 
wavy line in the diagram for 7/j corresponds to the multiphonon propagator 
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Pekarian distribution, and the vertex function transforms into a sum of phonon satellites 
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lK = ye-'—^H{nn). (5) 



Here S = v/Ml is the Huang-Rhys factor and v = \^ /Ml is the polaron shift. This 
equation is valid at kT < Ml and v > pt"^, otherwise the satellites smear into a single hump 
around the maximum of Pekarian distribution. Since HQ ^ E^, only one phonon replica 
satisfying the condition \A — h — nhQ\ ~ E^ survive at given magnetic field B . This 
means that changing B one may "scan" the Pekarian function in a certain interval. In case 
illustrated in Fig. [3] second system of sateUites arises with Kondo peaks satisfying condition 
I A + /i - (n + m)hn\ ~ Ek- 
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Coupling constant a = Ju/ep 



FIG. 5: Phase diagram illustrating the compe- 
tition between Kondo temperature and relax- 
ation damping. Solid line stands for character- 
istic Kondo energy. Dashed line corresponds to 
the damping parameter. The phonon- assisted 
Kondo tunneling is effective for a > Oc, where 
Uc ~ 10~^ for TMOC coupled with metallic 
leads. G/Gq is ehnancement factor for ZBA in 
tunnel conductance. 



Since phonon-assisted Kondo tunneling is essentially non-equilibrium process, one should 
estimate the contribution of decoherence and dephasing effects. Similar problem was dis- 
cussed in [qI, llQl] for a situation where the gap A is compensated by finite source-drain 
bias. To evaluate these effects, one should calculate the damping of S-T excitation (imag- 
inary part of the corresponding self energy stemming from the vertex part shown in Fig. 
S]). Performing calculations in analogy with p], one gets for the lifetime r the estimate 
U/t ~ {pJ^i^Q. This damping should be compared with the Kondo energy extracted from 
(jl]), Ek ~ ei?exp(— 1/pJ/j). The competition between two quantities reminds that between 
the indirect exchange /j„ ~ pJj^ and the Kondo energy in the Doniach diagram for Kondo 
lattices [7|. However, unlike the Doniach dichotomy, in our case Ex dominates in the most 
part of the phase space because h/r contains small parameter Q/D ^ 1 comparing to /j„ 
(see Fig. [5]) Thus, the spin-phonon relaxation is not detrimental for the phonon-assisted 
Kondo effect for realistic model parameters. One more interesting situation is the case 
where two vibration modes exist in a molecule, one of which may be put in a resonance with 
the state Ti at some field h = hi, whereas another one may be tuned to the state Tj at 
h = h2> h. (Fig. iD. 

To conclude, we have demonstrated that the multiphonon emission/absorption processes 
may initiate Kondo effect in TMOC with the ground singlet state in a situation where the 
electron tunneling from metallic reservoir to a 3d orbit of TM is accompanied by polaronic 
effect in molecular vibration subsystem. The release of vibrational energy compensates 
the singlet-triplet gap in the spin excitation spectrum. Varying the Zeeman splitting in 
external magnetic field, one may scan the Pekarian distribution of Kondo-phonon satellites 




FIG. 6: Two-mode regime. 

in tunneling conductance. Magnetic field tuning changes drastically the character of phonon- 
assisted Kondo screening. In zero magnetic field the spin S" = 1 is underscreened by phonon- 
assisted Kondo processes [3|], whereas the Zeeman splitting results in reduction of effective 
spin from 1 to 1/2 and hence to complete Kondo screening. One make hope to observe 
this phenomenon in the device with suspended nanotube, which plays part of a quantum 
dot HI]. The frequency of the stretching mode in this nanotube is comparable with the 
values of magnetic field used in standard experiments (^^strctch ~ TOO/ieV), so one may 
hope that the zero bias anomalies in Kondo tunneling assisted with one and two phonons 
(n = 1 in Fig. [3]) may be detected in this device. 
Authors are greatly indebted to M R Wegewijs and A. Hiittel for stimulating discussions. 
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